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Abstract. We present an algorithm for computing the integral closure 
of a reduced ring that is finitely generated over a finite field. 



Leonard and Pellikaan [1] devised an algorithm for computing the integral 
closure of weighted rings that are finitely generated over finite fields. Pre- 
vious algorithms proceed by building successively larger rings between the 
original ring and its integral closure, [21 [6l El [9j [Til Q2|- The Leonard- 
Pellikaan algorithm instead starts with the first approximation being a 
finitely generated module that contains the integral closure, and succes- 
sive steps produce submodules containing the integral closure. The weights 
in [1] impose strong restrictions; these weights play a crucial role in all steps 
of their algorithm. We present a modification of the Leonard-Pellikaan al- 
gorithm which works in much greater generality: it computes the integral 
closure of a reduced ring that is finitely generated over a finite field. 

We discuss an implementation of the algorithm in Macaulay 2, and pro- 
vide comparisons with de Jong's algorithm [2]. 

1. The algorithm 

Our main result is the following theorem; see Remark II .51 for an algorith- 
mic construction of an element D as below when R is a domain, and for 
techniques for dealing with the more general case of reduced rings. 

Theorem 1.1. Let R be a reduced ring that is finitely generated over a 
computable field of characteristic p > 0. Set R to be the integral closure 
of R in its total ring of fractions. Suppose D is a nonzerodivisor in the 
conductor ideal of R, i.e., D is a nonzerodivisor with DR C R. 

(1) Set Vq = jjR, and inductively define 

V e +i = {/ S 14 | f G V e } fore^O. 

Then the modules V e are algorithmically constructible. 

(2) The descending chain 

V DV 1 ^V 2 2V 3 D ■■■ 
stabilizes. IfV e = V e +i, then V e equals R. 
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The prime characteristic enables us to use the Frobenius or p-th power 
map; this is what makes the modules V e algorithmically constructible. 

Remark 1.2. For each integer e ^ 0, the module DV e is an ideal of R; we 
set U e = DV e and use this notation in the proof of Theorem 11.11 as well as 
in the Macaulay 2 code in the following section. The inductive definition of 
V e translates to Uq = R and 

U e+1 = {reU e \r p e D^Ue] for e > . 

Proof of Theorem \l.l\ (1) By Remark \1.2\ it suffices to establish that the 
ideals U e are algorithmically constructible. This follows inductively since 

U e+1 = U e n ker (R — R — 2— > R/D^Ue) for e > , 

where -F is the Frobenius endomorphism of R, and tt the canonical surjection. 

(2) By construction, one has V e +\ C V e for each e. Moreover, it is a 
straightforward verification that 

V e = {f eV \ f pt G y for each t < e} . 

Suppose f £ R. Then / pl G i? for each i ^ 0, so D/ pl G i?. It follows that 
/£ V e for each e. 

If = V'e for some positive integer e, then it follows from the inductive 
definition that V e +i = V e for each % ^ 1. 

Let ui, . . . , v s : R — * Z U {cxd} be the Rees valuations of the ideal DR, 
i.e., Vi are valuations such that for each n G N, the integral closure of the 
ideal D n R equals 

{r G R | Uj(r) ^ nvi(D) for each i} . 

Let e be an integer such that p e > Vi{D) for each i. Suppose r/D G V e . 
Then (r/D)P e G Vb, so rP e G D^~ X R. It follows that 

P e ^(r) ^ (p e - l)vi(Z)) 

for each i, and hence that 

Vi {r) > Vi{D)-Vi{D)/p e > Vi (D)-l 

for each i. Since Vi(r) is an integer, it follows that Vi{r) ^ Vi(D) for each i, 
and therefore r G Di?. But then r belongs to the integral closure of the 
ideal DR in R. Since principal ideals are integrally closed in R, it follows 
that r G DR, whence r/D eR. □ 

Remark 1.3. We claim that if R is an integral domain satisfying the Serre 
condition S2, then each module V e is S2 as well. 

Proceed by induction on e. Without loss of generality, assume R is local. 
Let x, y be part of a system of parameters for R. Suppose yv G xV e+ i for 
an element v G V e+ \. Then yv/x G V^+i, i.e., yu/x G V e and y p v p /x p G V^, 
or equivalently, yt> G xF e and y p v p G x p y e . Since V e is S2 by the inductive 
hypothesis, it follows that v G xV e and u p G x p V e , hence v G xV e +\. 
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Remark 1.4. In the notation of Theorem II .1\ suppose e is an integer such 
that V e = V e +i. We claim that the integral closure of a principal ideal aR is 

{r G R | D/ G a pl R for each i ^ e + 1} . 

To see this, suppose r is an element of the ideal displayed above. Then 
Dr p = gaP for some g £ R. Since 

D(r/af G R for each i < e + 1 , 

it follows that 

D{g/D) pl G R for each i < e . 
But then g/D G 14, which implies that g/D G V{ for each i. Hence 
D(r/a) pl G R for each i, equivalently r G aJ2. 

Remark 1.5. Let R be a reduced ring that is finitely generated over a 
perfect field K of prime characteristic p. We describe how to algorithmically 
obtain a nonzerodivisor D in the conductor ideal of R. 

Case 1. Suppose R is an integral domain. Consider a presentation of R 
over K, say = iT[zi, • • • , %n]/(fi, ■ • • , fm)- Set h = height^, f m ). 
Then the determinant of each h x h submatrix of the Jacobian matrix 
(dfi/dxj) multiplies R into R; this may be concluded from the Lipman- 
Sathaye Theorem ([5] or \10\ Theorem 12.3.10]) as discussed in the following 
paragraph. At least one such determinant has nonzero image in R, and can 
be chosen as the element D in Theorem I l.li Other approaches to obtaining 
an element D are via the proof of [101 Theorem 3.1.3], or equivalently, via 
the results from Stichtenoth's book |8j. 

Let J be the ideal of R generated by the images of the h x h submatrices 
of (dfi/dxj). We claim that J is contained in the conductor of R. By 
passing to the algebraic closure, assume K is algebraically closed. After 
a linear change of coordinates, assume that the Xi are in general position, 
specifically, that for any n—h element subset A of {x\, . . . , x n }, the extension 
K[A] C R is a finite integral extension, equivalently that K[A] is a Noether 
normalization of R. By the Lipman-Sathaye Theorem, the relative Jacobian 
Jr/k[A] is contained in the conductor ideal. The claim now follows since, as 
A varies, the relative Jacobian ideals Jr/k[A] generate the ideal J. 

Case 2. In the case where R is a reduced equidimensional ring, one may 
proceed as above and choose D to be the determinant of an h x h submatrix 
of (dfi/dxj), and then test to see whether D is a nonzerodivisor. If it turns 
out that D is a zerodivisor, set 

h = (0 : R D) and I 2 = (0 : R I x ) . 

Then each of R/I\ and R/I2 is a reduced equidimensional ring, with fewer 
minimal primes than R, and 

R = R/h x R/h ■ 

Hence R may be computed by computing the integral closure of each R/L L . 
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Case 3. If R is a reduced ring that is not necessarily equidimensional, 
one may compute the minimal primes Pi , . . . , P n of R using an algorithm for 
primary decomposition — admittedly an expensive step — and then compute 
R using Case 1 and the fact that 

R = R/P x x ••• x R/P n . 



2. Implementation and examples 

Here is our code in Macaulay 2 [3] , which uses this algorithm to compute 
the integral closure. 

Input: An integral domain R that is finitely generated over a finite field, 
and, optionally, a nonzero element D of the conductor ideal of R. 

Output: A set of generators for R as a module over R. 

Macaulay 2 function: 

icFracP = method(Options=>{conductorElement => null}) 
icFracP Ring := List => o -> (R) -> ( 

P := ideal presentation R; 

c := codim P; 

S := ring P; 

if o . conductorElement === null then ( 
J := promote (jacobian P,R); 
n := 1; 

detl := ideal(0_R) ; 

while detl == ideal (0_R) do ( 

detl = minors (c , J) ; 

n = n+1 

); 

D := detl_0; 
) else D = o . conductorElement ; 
p := char(R) ; 
K := ideal (1_R) ; 
U := ideal (0_R) ; 

F := apply (generators R, i-> i~p) ; 
while (U != K) do ( 
U = K; 

L := U*ideal(D~(p-l)) ; 

f := map(R/L,R,F) ; 

K = intersect (kernel f, U) ; 

); 

U = mingens U; 

if numColumns U == then {1_R} 

else apply (numColumns U, i-> U_(0,i)/D) 

) 

Since the Leonard-Pellikaan algorithm uses the Probenius endomorphism, 
it is less efficient when the characteristic of the ring is a large prime. In the 
examples that follow, the computations are performed on a MacBook Pro 
computer with a 2 GHz Intel Core Duo processor; the time units are seconds. 
The comparisons are with de Jong's algorithm [2] as implemented in the 
program ICfractions in Macaulay 2, version 1.1. 
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Example 2.1. Let F2[x, y, t] be a polynomial ring over the field F2, and set 

R = F2I/E, y, x 2 t, y 2 t). Then R has a presentation 

¥ 2 [x,y,u,v}/(x 2 v - y 2 u) , 

which shows, in particular, that x 2 is an element of the conductor ideal. 
Setting D = x 2 , the algorithm above computes that the integral closure of 
R is generated, as an i?-module, by the elements 1 and xyt. Tracing the 
algorithm, one sees that Vq is not equal to Vi, that V\ is not equal to V2, 
and that V 2 = V3. Indeed, these i2-modules are 

V = \r , Vi = -R + ytR , V e = R + ytR for e > 2 . 

X X 

As is to be expected, the algorithm is less efficient as the characteristic of 
the ground field increases: 

Table 1. Integral closure of ¥ p [x, y, u, v]/ (x 2 v — y 2 u) 



characteristic p 


2 


3 


5 


7 


11 


13 


17 


37 


97 


icFracP 


0.04 


0.03 


0.04 


0.04 


0.04 


0.05 


0.05 


0.13 


0.59 


icFractions 


0.08 


0.09 


0.09 


0.09 


0.14 


0.15 


0.15 


0.15 


0.15 



We remark that R is an affine semigroup ring, so its integral closure may 
also be computed using the program normaliz of Bruns and Koch pQ. 

Example 2.2. Consider the hypersurface 

R = ¥ p [u, v, x, y, z]/(u 2 x 4 + uvy 4 + v 2 z 4 ) . 

It is readily verified that R is a domain, and that t = ux^/v is integral 
over R. The ring R[t] has a presentation 

¥ p [u,v,x,y,z,t]/I , 

where I is the ideal generated by the 2x2 minors of the matrix 

fu t -z A \ 
yv x 4 t + y 4 J 

Since the entries of the matrix form a regular sequence in ¥ p [u, v, x, y, z, t], 
the ring R[t] is Cohen-Macaulay. Moreover, if p 7^ 2, then the singular locus 
of R[t] is V(t, y, xz, vz, ux) which has codimension 2, so R[t] is normal. 

If p = 2 then the ring R{t] is not normal; indeed, in this case, the integral 
closure of R is generated, as an i?-module, by the elements 

— ux + z^/mi vz + Xyfuv uxz + z 2 yfuv 

1, Vuv, , , . 

y y uy 

For small values of p, these computations may be verified on Macaulay 2 
using either algorithm; some computations times are recorded next. Here, 
and in the next example, * denotes that the computation did not terminate 
within six hours. 
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Table 2. Integral closure of ¥ p [u, v, x, y, z]/(u 2 x 4 + uvy 4 + v 2 z 4 ) 



characteristic p 


2 


3 


5 


7 


11 


icFracP 


0.07 


0.22 


9.67 


143 


12543 


icFractions 


1.16 


* 


* 


* 


* 



Example 2.3. Consider the hypersurface 

R = ¥ p [u, v, x, y, z]/(u 2 x p + 2uvy p + v 2 z p ) , 

where p is an odd prime. We shall see that R has p + 1 generators as an 
i?-module, but first some comparisons: 

Table 3. Integral closure of ¥ p [u, v, x, y, z\j (u 2 x p + 2uvy p + v 2 z p ) 



characteristic p 


3 


5 


7 


11 


13 


17 


19 


23 


icFracP 


0.07 


0.09 


0.27 


1.81 


4.89 


26 


56 


225 


icFractions 


1.49 


75.00 


4009 


* 


* 


* 


* 


* 



We claim that R is generated, as an i?-module, by the elements 
~z , and rf/PyiP-^/P for 1 ^ i < p - 1 . 



(2.3.i) i, \/r 

It is immediate that these elements are integral over R; to see that they 
belong to the fraction field of R, note that 

uy p + vz p 



V? 



xz 



' u(y 2 - xz)^- 1 )/ 2 
and that, by the quadratic formula, one also has 

'u\ 1 /p 



(2.3.2) 



-y ± vV 



xz 



Moreover, using (|2.3.2p . it follows that 

±{xu 1 'P + yv l /P) 



v x l p ^y~ 2 



xz 



and hence the i?-module generated by the elements (|2.3.ip is indeed an 
-R-algebra. It remains to verify that the ring 



A = R[y/y 2 -xz, u i/p v {p ~ l)/p 
is normal. For this, it suffices to verify that 



1 < i < p - 1 



B = R[y/ 



y 



xz, u 1/p V /p l 



is normal, since A is a direct summand of B as an A-module: use the grading 
on B where degx = degy = degz = and degu 1 ^ = 1 = degf 1 ^, in which 
case A is the p-th Veronese subring © iGN B{ p . The ring B has a presentation 
¥ p [x, y, z, d, s, t]/I, where / is generated by the 2x2 minors of the matrix 

r y + & z s 
x y — d —t, 
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and s i — ► u 1 ^, t t— > u 1 ^, d i— ► \J y 2 — xz. But then — after a change of 
variables — B is a determinantal ring, and hence normal. 
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for his feedback, and to Amelia Taylor for valuable discussions and help with 
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